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g;; : Abstract 

We investigate the spacetime of a thick gravitating domain wall for a 

general potential V{^). Using general analytical arguments we show that all 

^O \ nontrivial solutions fall into two categories: those interpretable as an isolated 

-Y-v \ domain wall with a cosmological event horizon, and those which are pure false 

^ ' vacuum de Sitter solutions. Although this latter solution is always unstable 



OO 



to the field rolling coherently to its true vacuum, we show that there is an 

(^ , additional instability to wall formation if the scalar field does not couple too 

CS \ strongly to gravity. Using the A<^^ and sine-Gordon models as illustrative 

Q"^ \ examples, we investigate the phase space of the gravitating domain wall in 

^^ ' detail numerically, following the solutions from weak to strong gravity. We 

O ■ find excellent agreement with the analytic work. Then, we analyse the domain 

I . wall in the presence of a cosmological constant finding again the two kinds of 

Vh ■ 
{^JQ, solutions, wall and de Sitter, even in the presence of a negative cosmological 

constant. 
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I. INTRODUCTION 

The study of topological defects has wide applicability in many areas of physics. In 
the cosmological arena, defects have been put forward as a possible mechanism for structure 
formation [jl|, and while recent work on global defects 0,|^ indicates that these were probably 
not responsible for structure formation, the intriguing discovery of a non-gaussian signature 
in the microwave background [^ leaves open the possibility that defects were around at some 
point in the evolution of the universe. 

Of all the topological defects, domain walls are the most deceptively simple to study. 
They correspond to solitons in 1+1-dimensions, which are extended in two spatial directions 
to form a wall structure. Because they depend on only one spatial coordinate, the distance 
from the wall, the solutions in the absence of gravity can often be written in closed ana- 
lytic form, and for certain potentials the models are completely integrable and the defects 
have the full interpretation of solitary travelling waves. In the presence of gravity however, 
the situation changes, gravity destroys the integrability of the theory, and also, except in a 
perturbative sense, the analytic nature of the solutions. Fortunately, since the walls corre- 
spond approximately to hypersurfaces in spacetime, there is a well defined way of analysing 
their gravity using Israel's thin wall formalism ^j, in which the wall is approximated by an 
infinitesimally thin hypersurface, and most of the literature concerning domain wall gravity 
uses this method. 

The domain wall is a rather interesting object gravitationally: unlike almost all of the 
other topological defects (the exception being the global string 0), its metric is not static 
but time-dependent [0,||, having a de Sitter-like expansion in the plane of the wall. Ob- 
servers experience a repulsion from the domain wall, and there is an horizon at finite proper 
distance from the defect's core. This horizon can be interpreted as a facet of the choice of 
coordinates, which usually use the flat space wall solution as a starting point, and impose 
planar symmetry on the domain wall spacetime. However, it is possible to use a different set 
of coordinates ^J^ in which the wall has the appearance of a bubble which contracts in from 
infinite radius to some minimum radius, then re-expands, undergoing uniform acceleration 
from the origin. The 'horizon' is then simply the lightcone of the origin in these coordinates, 
and is somewhat similar to the horizon of Rindler spacetime. 

The crucial physical difference of a domain wall spacetime, as opposed to that of the local 
cosmic string or monopole, is the presence of this cosmological horizon, which introduces 
a second length scale into the system. Ordinarily, a defect possesses one length scale, w, 
which is a measure of its thickness. However, the distance to the event horizon of the domain 
wall gives another length scale, z^, which can be compared to w. Since these lengths are 
given in terms of the coupling constants of the theory, taking a thin wall limit turns out to 
be a very artificial construction in terms of these underlying parameters, and the issue of 
the self-gravity of thick walls becomes more pertinent. After the original work by Vilenkin, 
Ipser and Sikivie 0J^, attempts focussed on trying to find a perturbative expansion in the 
wall thickness |]TD|,|TT| both for the purpose of discovering the motion of the wall, as well as 
verifying that the hypersurface formalism was a good approximation to the true gravitational 
field. These results backed up the thin wall approximation, the main difference being the 
presence of a sub-dominant tension transverse to the wall. 

With the suggestion of Hill, Schramm and Fry [|l^ of a late time phase transition with 



thick domain walls, there was some effort at finding exact thick solutions |T3|JT^, see also 



15| , [T6[| . However, such defects were supposed to be thick by virtue of the low temperature 
of the phase transition responsible for defect formation. The suggestion that Planck scale 
topological defects could be responsible for inflation [|l^|l8l then reopened the issue of thick 



domain walls. 'Thick' then means thick compared to their natural de Sitter horizon, and 
therefore this appears to be a very strong gravity situation. The issue of whether a domain 
wall can survive in a Friedmann-Robertson- Walker (FRW) universe with horizon size com- 
parable to wall thickness has been analysed in some detail for the case of no gravitational 
back reaction of the wall on the FRW background |n,5D|, and in the case of Euclidean in- 



stantons on a de Sitter background including self-gravity ||2T|J2^; however, to our knowledge, 
a systematic analysis of the strong self-gravity of thick domain walls with an arbitrary field 
theory potential has not been carried out. 

In this paper we perform a detailed analysis of the strong gravity of thick domain walls. 
With a combination of analytical and numerical results, we map out the parameter space 
in which a wall-like solution exists. Although we focus on two main examples of the sine- 
Gordon and A$^ wall, we also illustrate to what extent, and for what potentials, the results 
will also hold in general. We begin in the next section by setting out the general formalism 
for a thick domain wall, deriving the metric and Einstein equations, as well as defining what 
we mean by a 'wall' solution. We also show how to generalise a coordinate transformation 
which takes a planar domain wall spacetime with an horizon into a flat spacetime with an 
accelerating wall to the case of walls with finite thickness. The solution outside the wall 
horizon is shown to be related to an FRW cosmology with a slowly rolling scalar field. 

In the following section we derive analytic results for the thick wall. It turns out that 
there are two possibilities for the scalar field from which the 'wall' is constituted: it can 
either be a planar wall or a de Sitter solution. Since the de Sitter solution is always possible, 
but not necessarily stable, we derive analytic bounds on the gravitational coupling strength 
of the wall for when the de Sitter solution is the only possible solution, and when it becomes 
unstable to decay into a wall-like solution. Note that the de Sitter solution always contains an 
instability to the field rolling coherently down the potential well — we will not be interested 
in these instabilities, only in those which would lead to the formation of a domain wall. 
We derive these bounds for both our chosen field theory models, as well as for a general 
potential. We then provide, in the case of weak gravitational coupling, some perturbative 
solutions for the wall and its gravitational fields. These can be readily compared to existing 
results in the literature. We then present numerical results backing up and extending the 
analytic work. In the penultimate section we consider the domain wall in the presence of a 
cosmological constant, and conclude in the final section. 

II. PLANE-SYMMETRIC SPACETIMES 

We start by setting up the general framework for a domain wall coupled to gravity. We 
initially consider a general matter Lagrangian 

CM = {Va^y-U{^), (1) 

where $ is a real Higgs field and the symmetry breaking potential f/($), has a discrete 
set of degenerate minima. We assume that the spacing of these minima is proportional to 



the (dimensionful) parameter 77, which sets the symmetry breaking scale, and that U{^) is 
characterised by a scale Vp = U{^y), where $f is a local false vacuum situated between 
successive minima ($f = is a conventional choice). For example, in the usual 'kink' model, 
[/($) = A (<l>^ — Tj'^) , and we see rj directly, with Vf = Ar/^. For convenience, we scale out 
dimensionful parameters via 

X = $/r/, e = SnGr]^ (2) 

The dimensionless parameter e, which we call gravitational strength parameter, characterises 
the gravitational interaction of the Higgs field. Then, defining V{X) = U{riX)/VF,^ 



SnGCM = -^ 



w 



w^iVaXf-ViX) , (3) 



where w = ^/ ^^Qy represents the inverse mass of the scalar after symmetry breaking, and 

of course will also characterise the width of the wall defect within the theory. The equations 
of motion following from @ are simply 

Without loss of generality we can set w = 1 (which amounts to choosing 'wall' units rather 
than Planck units) and, looking for a static solution in fiat space, we see that (^ can be 
integrated directly to give 



X" = V{X), (5) 



which has an implicit solution 



/ = z- zo, (6) 

^F^/nx) 

where Xp = X{zo) is the false vacuum. For example, in the A^'* model above, z — zq = 
/q dX/{l — X^) = tanh~^X, and we get the usual kink solution centered on zq: X = 
tanh(2; — zq). Another model that we will be exploring in detail is the sine-Gordon model, 
with V'(<l>) = Ff[1 + cos($/?7)]/2, in which case z - zq = 21ntan(7r/4 + X/4). 

We now look for a plane-symmetric gravitating domain wall solution, since this represents 
the most obvious intuitive generalization of the fiat space domain wall. We will consider 
coordinate transformations of this solution at the end of the section. The metric therefore 
will have planar symmetry (i.e. Killing vectors dx, dy, xdy — ydx)-, and in addition will display 
reflection symmetry around a surface, 2; = say, which represents the location of the wall 
(defined by X = X-p)- If we choose z to be the proper distance from the wall, then the 
metric may be written as 

ds^ = A^{z)dt^ - B^{z, t) (dx^ + dy"^) - dz"^, (7) 



^Note that V{Xf) = 1 by definition. 



with the associated Einstein equations derived from the Lagrangian (|l]) coupled to gravity 
through the usual Hilbert term, Cq = —R/2e, being 

Rab = 2eX,aX,,-eViX)gat. (8) 

Before writing these equations explicitly, we will first examine what we mean by a domain 
wall solution, since this will require some rather specific boundary conditions at z = 0, and 
for large z. We will define a wall solution to be a function X{z) of the proper distance 
from the wall which at z = is at a local maximum of V{X), Xp, and which falls towards 
distinct minima on either side of the wall. Assuming that V{X) is locally symmetric around 
the maximum, then {X{z) — Xp) will be an odd function of z. This restriction embodies 
the idea of a defect, in which the field falls to distinct vacua on either side of its core, and 
settles to a topological, rather than radiative, configuration. It is possible that the only 
nonsingular solution satisfying these criteria is X = Xp, in which case the spacetime will be 
de Sitter; indeed, this is always a possible solution to the equations of motion satisfying the 
above criteria, though it will not necessarily be stable. Clearly, there is always an instability 
corresponding to a coherent roll of the field towards one of its vacuum values, however, since 
we are interested in spacetimes with the interpretation of a domain wall, we will ignore this 
mode, and say that de Sitter spacetime is 'stable' if there is no odd perturbation of {X — Xp) 
which is growing in time. As we will see, the scalar field does not fall all the way to its 
minimum value within the range of validity of the coordinates in (|^), which is why we do 
not place any specific boundary conditions on X for large z, however, it will turn out that 
we require X'{zfi) = for a nonsingular solution. Finally, we can choose t to set A{0) = 1, 
and reflection symmetry requires ^'(0) = 0. 

Turning to the Einstein equations (§), we see that 

BA' B' 

which implies B = b{t)A{z). Then the relation 

Rl-R: = A^[y^-l^=0 (10) 

yields 

bit) = e'\ (11) 

so that the equations of motion for the gravitating wall finally reduce simply to 




3 

^ 1 

A2 + 3 



2X'2 + V{X) (12b) 



X'^-V{X) . (12c) 



Note that since A" < 0, once A' becomes negative it will always be bounded away from 
zero, therefore there is some finite Zh for which A(zh) = 0, and we have either a physical or 



coordinate singularity. Thus we see immediately that there is no nonsingular solution with 
k'^ = 0. It is also easy to see that the false vacuum de Sitter solution is given hj X = Xp, 
A = cos kz, e = e/3.g 

Determining therefore whether or not a wall solution exists reduces to investigating the 
system of equations (|13). Clearly, for small e we might expect a wall solution to the above 
equations to exist, and to be given by a perturbative expansion around flat space. For large 
e, since A' = 0(e), the A'X'/A term in ( |12a| ) will drive the solution to a singularity for 
nonzero {X — Xp), hence for large e we expect only the de Sitter solution. For intermediate 
values of e, except for some very special cases, an analytic solution does not exist and we 
have to numerically integrate the equations. 

Before proceeding with the details of this analysis, we conclude this section by comment- 
ing on a coordinate transformation of the plane symmetric metric (^ , which transforms the 
defect into an accelerating bubble wall. Defining 

X* = A{z)e'''x (13a) 

y* = A{z)e''y (13b) 

t* ~ z* = ~lA(z)e''' (13c) 

k 

t* + z* = ^A(z)e-''' - kix^ + y2)A(;z)e'=* (13d) 

k 

gives an alternate form of the line element: 

ds'^ = dt*^ - dx*^ - dy*^ - dz*^ +l^-l\dz\ (14) 

where z is given implicitly by 

k\t*^-^*^) = -A\z). (15) 

The wall (i.e. the zero of X — Xp) is located at x*^ — t*^ = l/Zc^, and as we will see, for small 
values of the gravitational coupling e. A' rapidly approaches k = 0(e) outside the core of 
the wall. Therefore, we see that the spacetime in these coordinates is approximately fiat, 
with the wall being located at a spacelike hyperboloid at a distance 1/k from the origin. 
This corresponds to a wall undergoing uniform acceleration, contracting in from infinity, 
reaching a minimum radius, and re-expanding outwards. The horizon (^(^^h) = in the 
old coordinates) is now the lightcone centered on the origin, and the region exterior to the 
horizon is the causal future and past of the origin. Thus the coordinate transformation (p!3[) 
generalises the thin wall transformation given in [^, and discussed in detail in section three 

ofi- 

For larger values of e the spacetime will not be fiat, and X will not be at its true 
vacuum value near the horizon. In particular, setting t* = Tcoship, x* = r sinh ■?/' sin ^ cos 0, 



^This somewhat less familiar form is merely one of the many coordinate transformations of de 
Sitter, and can be reduced to the more familiar form ds"^ = dr"^ — e^^'^d'x^ via the transformation 



e = e COS kz, C, = tan(A;z)e /k. 



etc., then gives a rather famihar form for the metric exterior to the horizon, i.e. inside the 
hghtcone of the origin: 



ds^ = C%T)dT^ 



T 



dip"^ + sinh^ ip{de'^ + sin^ 



(16) 



the metric of an open FRW universe. Thus if the scalar field is not at its vacuum value at 
the horizon, its evolution outside the horizon is given by the rolling of a scalar towards its 
minimum in an open FRW model, a well studied problem! 

III. SELF-GRAVITATING DOMAIN WALLS 



We now want to find solutions to ([T^) representing an isolated domain wall for various 
values of the parameter e. We start by considering a general potential, V{X), finding an 
analytic perturbative solution for small e, and showing that if e is sufficiently large there is 
no wall solution. We also demonstrate an instability of the false vacuum de Sitter solution to 
wall formation for e < emax, where emax depends on the second derivative of the potential at 
the origin. We then explore the particular cases of the A^^ kink and the sine-Gordon model 
in some detail, first making reference to the analytic work, then describing the solutions 
for intermediate values of e with the help of numerical work. Both analytic and numerical 
work show the presence of a phase transition in the behaviour of the solutions between wall 
existence and nonexistence. 

It is clear that when the gravitational strength parameter is set to zero, one gets the flat 
space solution A = 1, with X being given by (^). Let us now consider small values of e, 
typically e <^ 1. Then we can expand the fields X and A in powers of e 



X = Xo + eXi + (e') 
A = Ao + eAi + (e^ 



(17a) 
(17b) 



where X„, An are now independent of e and Ao,Xo are the fiat space solutions. 
The field equations (|12D to lowest order in e give 



A'l 



2X',' + ViXo 



Vl2 
-^0 



1 9V 



^o(^) 



The boundary conditions of A and X give 

Ai{0) = A[{0) = 0, Xi(0) = 0, Xi ^ for large z; 
( |18a|) and ( |18b|) can be integrated to give 

dX 



A, 



V{Xo)dz 



WdX 



(18a) 
(18b) 

(19) 

(20) 



and 




Xi = -7;XU^2 A?--^]dz, (21) 



which can also be expressed as an iniphcit integral in terms of V{X) and X. Finally, noting 
that 

Xl,X[ - X,X'^ = ^Af + X^Xm, (22) 



^ |12c| ) implies 



2 



3 



2e 



2 



A? + -(XiXj - X^XO = -— X^(0)X((0). (23) 



3 



Before exploring these solutions for specific models for ^(X), we can make some general 
statements about the self-gravitating wall. First of all, since A" is strictly negative, A will 
asymptote 1 — kz, where k is 0(e) from the above relation. This means that Ai will cease 
to be small at a distance of order e~^ from the wall, and our expansion procedure strictly 
speaking breaks down. However, since A[ is not growing, it is clear that continuing the 
expansion to higher orders in e will merely produce minor corrections to A, and will not 
alter the qualitative behaviour, namely, that A{z) has a zero at a distance of order e~^ from 
the wall (see figure 0). Since gu = A'^ this is simply the event horizon which is familiar from 
the thin wall approximation. 

As e increases, the effect of the scalar field's energy-momentum on the geometry increases, 
and we expect the horizon to move closer to the wall, roughly until e = 0(1). For large e, the 
expected horizon would be well inside the core of the wall and two possibilities now emerge. 
First of all, the scalar field could simply ignore the geometry, and fall away minutely from its 
false vacuum, causing an horizon at some small value of z. Since the spatial gradient of such 
a solution would be relatively small, the energy-momentum would be vacuum dominated, 
and we would expect the spacetime to be very close to the de Sitter solution. Outside 
the horizon, the scalar field would roll to its vacuum value as described in the last section. 
The alternative is that there is a phase transition in the behaviour of X, that is, that X 
either has some nontrivial odd form, approaching reasonably close to its vacuum value at 
the horizon, or X = Xp. In other words, X must either roll significantly away from Xp, 
or not at all. There are two reasons why we expect the latter scenario to hold. The first is 



that Basu and Vilenkin [^ observed just such a phase transition in studying the problem of 
wall instantons. The other reason for suspecting a phase transition lies in the behaviour of 
field theory solutions on compact surfaces. For example, two of the authors have studied a 



cosmic string interacting with the event horizon of an extremal black hole ||2^ . There, there 
are nontrivial solutions with the string piercing the horizon while the string fields can fall 
reasonably close to their vacuum values around the horizon; however, there is a transition 
when the string becomes sufficiently thick relative to the black hole: the event horizon can 
no longer support a nontrivial solution and the flux of the string is expelled — the only 
horizon solution is the trivial one. 

First of all, let us show that for small e the false vacuum de Sitter solution is unstable. 
Recall that the de Sitter solution is 

ds^ = cos\kz)dt^ - e^^* cos\kz){dx^ + dy^) - dz^ (24) 



with X = Xp and k"^ = e/3. This solution will be 'stable' (i.e., stable to wall formation) 
if there is no perturbation of X — Xp which is an odd function of z and growing in time. 
Setting X = Xp + ^, and noting that corrections to the geometry are 0(^^), we see that an 
instability must satisfy the time-dependent linearized perturbation equation 



e" - 3kt8in{kz)^' - sec\kz)[^ + 2k^] - -^ 



l.d^V 







Xv 



(25) 



2^dX^ 

(where a dot indicates differentiation with respect to t). This equation does indeed have 
unstable solutions for /c^ = e/3 < |V""(Xf)|/8, the dominant instability being given by 

^ = e''''*sin/c2;(cos/c2;)' 

with 

1 



(26) 



_5 

~2 ^ 2 



Thus for e smaller than 



9 - 2V"{XF)/k^. 



\v"ix. 



(27) 



(28) 



the de Sitter solution is unstable to wall formation. 

Now let us examine whether a wall solution can exist for large e. For a nontrivial wall 
solution, we require X'(0) > 0, and for a nonsingular solution, we require X'{zi^) = 0. Now 
taking the derivative of the field equation ( |12aD we get 

' 1 d^v' 

+ 



X" 



-3^X" + X' 
A 



-3- 



A" 



A 




2 9X2 



A' 



X'F(z) 



where F[z) may be rewritten using (|T^ as 



Fiz) 



3eX 



/2 



k^ id^v 

_|_ '3. I 

A2 2 9X2 



Now, at 2 = dm) gives 3A;2 = e[l - X'(0)2], hence 



X'"(0) = X'(0) 



3e - Qk' 



\v'\Xp) 



> X'(0) 



e+2^"(^F) 



(29) 



(30) 



(31) 



Therefore, if e > |\^"(Xf)|/2, X'"(0) > 0, and X' is increasing away from z = 0. Moreover, 
if |\^"(X)| is maximized at Xp, then F{z) is strictly increasing away from z = and X' 
can never be zero at the horizon. Thus with a minimal assumption on the nature of the 
potential, we have shown that there is no wall solution possible if 



e > - \V"{Xp) 



(32) 



This argument does not make any assumptions as to the behaviour of the geometry, it simply 
relies on general properties of the potential. 

To reiterate, we have shown that for e < emax there are two solutions: a wall spacetime 
and a de Sitter solution, the latter being unstable to wall formation. For e > emax, the de 
Sitter solution is 'stable', and for e > 4einax/3, the domain wall solution can analytically 
be shown not to exist. To determine whether de Sitter is the only solution for e^ax < e < 
4emax/3, the problem must be examined numerically. We now do this, and obtain explicitly 
the perturbative solution for the A$^ and sine-Gordon potentials. 



A. The A$^ kink 



In this case the potential V{X) is given by V{X) = {X^ — 1)^, and in flat spacetime 
(0(6")) we have the usual flat domain wall solution plotted in figure |l|, 



Xq = tanh(2;). 



(33) 
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-2 - 

-3 F- 



FIG. 1. The flat spacetime solution. The solid line is X = tanh(z), the solution for the 
Goldstone model, and the broken line is X = 4 arctan ( e 2 j — vr, the solution for the sine-Gordon 
case of section [II B. 



Integrating (^), (0) and calculating 

k = 



one gets, 



z -\ — tanh z 
3 



-sech^z 
2 

2, , 1 ,2 
- log cosh z tanh z 

3 ^ 6 

|. + o(.^). 



(34a) 
(34b) 
(34c) 



which represent the first order gravitational corrections to X and A. Note also that Xi does 
indeed satisfy the boundary conditions (plQl). The distance to the event horizon is given by 
Zh — 3/2e. Putting together our results we get to order 0(e) 



X 
A 



tanhz sech z 

2 



z -\ — tanh z 
3 



1-'- 
3 



2 In cosh z -\ — tanh z 
2 



(35a) 
(35b) 



This solution is compared on figure ^ with the one found numerically, for e = 0.1. 

For larger values of e, we must resort to numerical methods to find solutions of (|12|). 
Here, we have used the routine solvde from |2^. The wall solutions that we obtain are 
qualitatively the same as the one shown on figure ^ Note that as mentioned previously, X 
does not go to its asymptotic value at z^ (and consequently that the energy density does not 
tend to zero at the horizon). In fact, we do not solve ([T^) as written in section |I|; instead, 
we rewrite equation (|12b[ ) as 
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FIG. 2. Comparison between the solution obtained numerically (solid line) for e = 0.1 and the 
series to order 0(e) [equation (|35D]. (The two solutions for X appear identical at this scale.) 




2X'^ + V{X) 



0. 



(36) 



The system ( |12a| , |12b|) can now be written as three coupled first order ordinary differential 
equations (ODE's), 



X' = Y 

Y' = -3YZ 



I dV 
2dX 



Z' = - 



2Y^ + V{X) 



-Z' 



(37a) 
(37b) 
(37c) 



where Z = A' /A. These equations were solved for the boundary conditions X(0) = Z{0) = 

F(^h) = 0. 

To determine whether there is indeed a phase transition in the behaviour of the solutions, 
we examine the evolution of the value of the Higgs field at the horizon, X\z^^, as a function 
of e. For a wall solution, this value represents the maximum of the function X (see for 
instance figure |^a), and a value X\z^^ = Xp = corresponds to the the false vacuum de 
Sitter solution. We expect therefore that X\z^^ will drop from 1 at e = to for some 
e in the range [cmax, fcmax]- According to our previous discussion, it is the solutions in 
an intermediate range of e which interest us most. We find (figure ^a) that the scalar 
field undergoes a phase transition at the value e = emax = 3/2, in perfect agreement with 
the prediction (^Sf ), at which point wall solutions cease to exist, and only the de Sitter 
configuration remains. 

Figure ^b shows the evolution of the proper distance to the horizon as a function of e; for 
small e this approaches the first order prediction z^ ~ 3/2e (dashed line), but higher order 
corrections rapidly spoil the agreement. The proper distance to the horizon at the phase 
transition can be predicted by the condition cos(fczh) = 0, which — with k^ = emax/3 = 
1/2 — implies Zh = T^I\P1 ^ 2.221. 

Thus we see that the numerical work confirms the general analytic derivations given 
earlier in the section, and indicates that at emax = 3/2, the domain wall solution disappears 
entirely. 
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FIG. 3. Numerical solution of the equations ( [l^ ) for the A<^^ model. This solution was obtained 
for e = 0.9 (in which case the horizon was situated at a proper distance Zh = 2.789). The figure 
shows (a) the Higgs field, (6) the energy momentum tensor T^^ = T^y = T^t and T^^, (c) the 
function A[z) and (d) the metric component gtt{z) = ^^(-2)- 
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FIG. 4. (a) The evolution of X\z^^ as a function of e. (&) Log-log plot of the proper distance 
to the horizon as a function of e (solid line) compared with the first order prediction of z^ = 3/2e 
(dashed line). The dash-dotted line indicates the phase transition at e = 3/2,X|^j^ « 2.221. 
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B. The Sine-Gordon Potential 

Consider now the periodic sine-Gordon potential, V{X) = |(1 + cosX) = cos^X/2. As 
before to zeroth order in e, one gets the usual sine-Gordon soliton 

Xo = 4 arctan (e^) - vr, (38) 

as shown in figure |I|. 

Making use of (pO| , pl|) , we obtain the gravitational back reaction to order 0(e), 

A = l-Ae In cosh - (39a) 

X = 4 arctan ( e 2 ) — n — 6e z sech- (39b) 

and k = 2e + (e^), which fixes the horizon distance to Zh — l/2e. Note the agreement with 



equation (3.19) of Widrow's paper |jTT|, who also considered a sine-Gordon domain wall, 
although he did not compute the correction to the scalar field. 

Again, we must turn to numerical methods to find solutions for higher values of e. The 
results we find are qualitatively very similar to those obtained for A$^. In particular we 
observe again the phase transition predicted in the previous section for a general potential. 
This time, however, the analytic results predict emax = 3/16 and z^ = 7i/2k = 27r; again, 
this is in excellent agreement with the numerical results. 

IV. DOMAIN WALLS WITH A COSMOLOGICAL CONSTANT 

In this section we consider the previous theories in a universe with a non-zero cosmological 
constant A, which would correspond to gravitating domain walls in an infiating universe. 
The effect of this constant can be readily taken into account by modifying equations (|T2p as 
follows: 

= -- [2X'2 + V{X)] - -A (40b) 

= ^ + ^ [X" - V{X)] - U. (40c) 

There are two qualitatively distinct cases: if A > 0, the wall is embedded in a de Sitter 
background, whereas if A < it is in an anti-de Sitter background.^ The latter is of 




^Strictly speaking, for an anti-de Sitter background, in order to have the reflection symmetry 
around z = we need to have k'^ < 0, which for a real metric would require b(t) = cos kt. This in 
turn requires the {x,y} sections to be hyperboUc (see for example ^5|); however, since this does 
not affect the equations of motion for A{z), we will not discuss it further, and instead refer the 
reader to [^] (and references therein) for a detailed review of anti-de Sitter domain walls. 
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particular interest because the effect of the cosniological constant should counteract the 
(effective) cosniological constant created by the wall's back reaction. 

First let us review how the analytic arguments of the previous section are affected by a 
cosniological constant. Note that a false vacuum solution X = Xp will now have an effective 
cosniological constant Aog = A+e, hence the metric will be of the form p^ ) with k'^ = Aeflf/3. 
Therefore, the previous arguments go through essentially unchanged, but with Aefr instead 
of e. Therefore 



\V"{X, 



A 



(41) 



Obviously, the range of the instability is increased for negative A, and decreased for positive 
A. If A > ||y'(XF)|, then the false vacuum solution is 'stable' (and indeed the only one). 

This is illustrated in figures |^ and ^ which show the evolution of X\z^^ and z^ with e 
for both A$^ and the sine-Gordon models as well as for several values of the cosniological 
constant, A = —0.3, —0.2, . . . , 0.3. In partiular, note that for sine-Gordon (figure ^6) the 
formula above tells us that for A > 3/16 = 0.1875 the only solution is X = Xp; this is why 
we do not see the curves for A = 0.2 and 0.3. Note as well that the value of Zh at which 
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FIG. 5. Evolution of X|^^ in function of e and A (from right to left, A = -0.3, -0.2, . . . 0.2, 0.3). 



(a) shows the A$ case, and (b) shows the sine-Gordon case. In (b) we have actually divided X\ 
by vr to help the comparison with case (a). 



Zh 



the phase transition occurs (7r/\/2 for A<l>^, and 2n for sine-Gordon) remains unaltered by 
the inclusion of the cosniological constant, as expected from the discussion above. In fact, 
all the de Sitter solutions remain identical if e and A are allowed to vary but Aeg remains 
constant. (This is obviously not the case for the wall solutions, as e then multiplies terms 
containing the Higgs field.) 

Now let us turn to the solution in the anti-de Sitter case, A < 0. We now find three 
qualitatively distinct solutions. For very small e, the wall's self-gravitation cannot compete 
with the anti-de Sitter expansion and A' /A is strictly positive; in fact, it is easy to check 
that the solution plotted on figure |^a is A{z) = cosh(-i/|A|/3z). As one increases e, the 
potential is observed to decrease close to the wall's core, whereas the Higgs profile is slightly 
smoothed (figure ^b). In fact, this is the beginning of a complete change in the metric 



14 




20 



15 - 



Zh 10 



- \ 


\ 

\ 
\ 


\ 

\ 
\ 
\ 

\ 


A - -0.3 

\ 

\ - 


^==s^ 


=^^^^== 


^^^^1::;^:^ 


A = +0.3 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
e 

(b) 



FIG. 6. Distance from the wall to the horizon, as a function of e and for the same values of A 
as in figure ^. (a) was obtained for A<1>^ and (b) for sine-Gordon. Again, the broken lines show the 
values of Zh at the phase transition. 



function A{z): as the wall's gravitational interaction is switched on, A assumes the shape 
of a "double well," with a local maximum at the imposed boundary value A{0) = 1 and 
two local minima symmetrically situated at A{±z-[a) for some z^- As e increases, this 
double well becomes deeper, whereas z^^ moves away from the wall. Notice that so far the 
function A{z) is strictly positive, and therefore none of these solutions exhibit an event 
horizon. Eventually, however, for some critical value, ec, of e, the two minima of A{z) 
vanish as 2;m — > cxd (figure 0c). This would appear to be a thick wall version of the type 



II extreme domain wall spacetime of Cvetic and Griffies p^l, and is therefore presumably 
supersymmetrizable. For e > ec, the metric becomes negative at a finite distance z^^, thus 



giving rise to the wall's horizon. The Goetz solution |]13[ lies in this range. 

Figure |^ shows the parameter space (A,e), and the different kinds of solution that we 
find. It is interesting to note that the two lines separating the three phases seem to run 
parallel to each other in both cases, indicating that a phenomenon similar to the triple point 
observed in the phase diagram of water never occurs. This is to be expected, since as long 
as the wall does not have an event horizon it is constrained to take its asymptotic value 
at infinity. Of course, this topological constraint does not imply that the lines are parallel, 
merely that they cannot meet in the physical range e > 0; figure |^ then shows that the range 
of the parameter e over which the value of the Higgs field at the horizon is allowed to drop 
from 1 to is fixed. 



V. DISCUSSION 

To summarize, we have investigated the spacetimes of thick, gravitating domain walls in 
detail, using both analytical arguments and numerical integration. Both methods demon- 
strate the existence of a phase transition in the nature of the 'wall' solution, from being 
wall-like to a pure false vacuum de Sitter spacetime. We find that for walls much thinner 
than their cosmological horizon, the Vilenkin thin wall solution is a good description of the 
spacetime. For thicker walls the spacetime differs more markedly until finally there is only 
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FIG. 7. Solutions ^(2;) (solid lines) and Ai^z) to the sine-Gordon equations for A = —0.3 and 
e = (a); 0.2 (6); 0.367. . . (c) and 0.4 (c?). 
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FIG. 8. Parameter space (e, A) and the types of solutions found. 
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the de Sitter solution. The transition occurs abruptly, when the gravitational strength of 
the wall is order of magnitude unity. For walls less strongly gravitating than this critical 
value, the spacetime has the appearance of a gravitating domain wall, possibly lying within 
a background de Sitter or anti-de Sitter universe. For very strongly gravitating domain 
walls however, the de Sitter universe is the only solution with the symmetries we have im- 
posed, the cosmological constant being provided by the false vacuum energy of the wall. Of 
course this does not prove that the de Sitter solution is the only solution to the coupled 
Einstein-scalar field equations; by transforming to a global coordinate system for de Sitter, 
one can find a complete set of solutions for the tachyonic wave equation, hence we can find 
an instability of the required parity which is given in the notation of section |ITI| by 



C = sm{kz)[v{x,y,z,t)]-^''+'>^/\F, 



u+ 1 u 



5 7 ■ 

— ,z/+ -\v 



2' 



where v is given by (|27|) , and v by 



v{x,y,z,t) 



1 + cos^ kz (sinh kt + \k'^{x^ + y'^)e''^y 



(42) 



(43) 



and F is a hypergeometric function. However, since this instability depends on the {x, y} 
coordinates, it does not correspond to decay to a solution of the type considered in this 
paper, and presumably corresponds to decay to a time dependent spherical wall of the type 
found numerically by Sakai et. al. |^ which are of relevance to topological inflation ||T7| , p!8 



These solutions are not solitons in the sense of having a fixed profile in time, and therefore 
are outside the consideration of this paper^ 

This phase transition behaviour of the scalar is reminiscent of the flux expulsion of a 
vortex core by the horizon of an extremal black hole |2^j2^. There, it was the fact that 
in the extremal limit the black hole horizon 7i decoupled from the exterior spacetime that 
allowed a partially analytic analysis of the vortex equations on the S"^ surface that was the 
event horizon. The findings in that case, |^, parallel our results here very closely. There 
was always a solution corresponding to the fields taking their false vacuum values on TC, 
but for sufficiently thin vortices (relative to the horizon radius) there was an additional 
solution corresponding to a vortex anti- vortex pair at opposite poles (in the case of a string 
threading the black hole). In that case however, the false vacuum solution for thin strings 
could be shown not to extend to a full solution in the exterior spacetime. Here, we always 
have a false vacuum de Sitter solution, which is unstable to wall formation, as well as the 
defect solution for low gravitational coupling. It is easy to see the common feature in these 
two problems — the compact nature of the spatial section upon which the defect must live. 
Defects on compact spaces have been analysed; for example. Avis and Isham [^ explored 
some years ago the A$^ solutions on a circle. They found exact solutions for the scalar field 
in terms of elliptic functions, and no solution other than false vacuum if the radius of the 
circle was too small. However, the crucial difference of our work to [^ and [^ is that 
we are not looking at defects on a fixed background, but looking for self-gravitating wall 
solutions without specifying their topology ab initio. 



We would like to thank Alex Vilenkin for discussions on this point. 
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The topology of a black hole event horizon is obviously compact, however, it turns out 
that in fact the topology of a domain wall spacetime is also compact [P]. To see that not 
only de Sitter spacetime but also the domain wall spacetime is topologically S^ x IR, consider 
the coordinate transformation ([l3D . If we define a fifth coordinate, w*, by 




w*= \l--^dz (44) 



then the wall metric becomes a slice of a five dimensional flat metric, and we can view our 
spacetime as a four- dimensional hypersurface embedded in five-dimensional flat spacetime 
in an analogous fashion to the de Sitter hyperboloid. From (0) the equation for this 
hypersurface is 

r2-x*2 = -^l!M. (45) 

For example, the de Sitter solution is A = cos kz, therefore w* = j: sin kz from (^^, and 
(^) reduces to t*^ — x*^ — w*"^ = — 1/fc^, the de Sitter hyperboloid. For small e on the other 
hand, the sine-Gordon wall from (|39D gives 



r^ z z 

w* = sech-dz = 4arctane2 — vr = Xq{z), (46) 

Jo 2 

hence the hypersurface is given by a hyperboloid which has been deformed by squashing in 
the w* direction, 

t*'-x*^ = -^(l + 4elogcos^) • (47) 

The spatial section is depicted in figure ^ which shows the t* = z* = slice to 0(e) for 
e = 1/30. 




FIG. 9. The spatial section of the weakly gravitating sine-Gordon domain walk The surface 
shown is the t* = z* = surface for e = 1/30. 

The spatial geometry of the domain wall is therefore topologically S^, and is similar to 
a discus, although the upper and lower surfaces are fiat almost to the edge of the discus. 
This corresponds exactly to our intuitive idea of the spacetime exterior to the wall being 
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FIG. 10. Bifurcation diagram for the gravitating domain wall. We plot here the normalized 
action 5 as a function of e. Solid lines indicate the stable solutions found numerically, and the 
dashed line represents the unstable de Sitter solution. 



flat, with a highly localised region of curvature generated by the wall itself located at the 
rim of the discus. It fleshes out the thin wall description of the spacetime, smoothing over 
the distributional singularity of the thin wall hypersurface. The two length scales of the 
domain wall spacetime correspond to the two radii of the discus, and as e increases, the 
discus radius shrinks, with its height remaining much the same, until the geometry is almost 
spherical, at which point the radius becomes too small to support a defect solution and 
we make the transition to a false vacuum de Sitter hyperboloid, with an exactly spherical 
spatial geometry. What is interesting in comparison with the Avis and Isham scenario is 
that it is the self gravity of the domain wall which produces a compactification of spacetime 
at its own characteristic scale, which ultimately becomes too small to support the defect 
itself. 



Finally, in analogy with |^ , we can make a plot of the normalized action for the domain 
wall solution versus the de Sitter solution 



S= f Cg + Cm = V^ [ V{X) e^'^'A^d^x = N f" A^V{X)dz 



(48) 



where A^ is a normalization factor, and there are no boundary terms from the gravitational 
part of the action. For the false vacuum de Sitter solution 5* = S/N = 2/3/c = 2/v^, and 
for the weakly gravitating A$^ model, S = 2/3. Figure |10| shows a plot of the action of the 



A$ wall against the false vacuum de Sitter solution, which indicates clearly the instability 
of the latter solution to wall formation for e < 3/2. (This can be compared to the instanton 
action plot obtained in pl| .) 
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